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Experiments and numerical simulations show that the energy spectrum of the magne- 
tohydrodynamic turbulence in the quasi-static limit deviates from Kolmogorov's k~ 5 ^ 3 
spectrum as the external magnetic field, or equivalently the interaction parameter, is 
£\j I increased. To explain this phenomena, we construct an analytical turbulence model 



with variable energy flux that arises due to the Lorentz-force induced dissipation. 



p*V The energy spectra computed using our model for various interaction parameters are 

in qualitative agreement with earlier experimental and numerical results. 
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I. INTRODUCTION 

Magnetohydrodynamic (MHD) flows are observed in two kinds of fluids: plasma and 
liquid metals. The flows in the stars and galaxies involve MHD plasma, while the interiors of 
the planet including the Earth, fission and fusion reactors (e.g., International Thermonuclear 
Experimental Reactor), industrial processes like metal plate rolling and crystallization etc. 
have liquid metals as the operational fluid. An important characteristic of the liquid metals 
is its low magnetic Prandtl number Pm, which is the ratio of the kinematic viscosity v 
and the magnetic diffusivity r\. Consequently, in many situations, especially in industrial 
applications, the magnetic Reynolds number Rm = UL/rj (U, L are the large scale velocity 
and length scales respectively) is quite small. Therefore, liquid metal flows under an idealized 
limit Rm — > 0, which is referred to as "quasi-static limit", is of major interest to scientists 
and technologists. 

In the quasi-static limit, the MHD equations get simplified significantly. Here, the induced 
magnetic field tends to be very small because of the large magnetic diffusivity, and it gets 
slaved to the velocity field that yields the Lorentz force as 

„ ° B l 1 d 2 " (1) 



p V 2 dz 2 

where p is the density of the fluid, u is the velocity field, and B = B z is the external 
uniform magnetic fielcU^. The Lorentz force induces additional dissipation of the kinetic 
energy. We will show in the present paper that the dissipation due to the Lorentz force 
reduces the inertial-range energy flux, and makes it scale-dependent. As a result, the energy 
spectrum gets steepened as the interaction parameter, which is the ratio of the Lorentz force 
and the nonlinear term pu • Vu, is increased. 

Several experimental and numerical simulations have been performed to study quasi- 
static MHD turbulence under a strong external magnetic field (see Knaepen and Moreau^ 
and references therein). Kolesnikov and Tsinober- and Alemany et al— performed exper- 
iments on mercury for low Rm. They reported that the energy spectrum for the velocity 
field is k~ 3 when the interaction parameter is significant. Kit and Tsinoner- argued that 
the observed spectrum is due to the two-dimensional nature of the flow, while Alemany et 
a/.- obtained k~ 3 spectrum by equating the nonlinear time scale with the Joule time scale, 
which was assumed to be wavenumber-independent. Sommeria and Moreau^ provided phe- 



nomenological arguments for identifying the time scales when the MHD turbulence with low 
magnetic Reynolds number becomes two-dimensional. 

Numerical simulations have been applied to study liquid metal flows under the quasi-static 
approximation or for low magnetic Reynolds number. Many numerical simulations^^— 
show steepening of the energy spectrum for large interaction parameter, similar to those 
seen in the experiments. It has been observed that for large interaction parameters, the flow 
becomes anisotropic with energy concentrated near the plane perpendicular to the external 
magnetic fielcr^^— . Ishida and Kaneda^ studied studied the anisotropic velocity correlation 
spectrum using theoretical and numerical tools. Favier et al.— argued that the flow is not 
entirely two-dimensional due to the presence of all the three components of the velocity field 
near the equator. This kind of fluid configuration is called "two-and-a-half-dimensional" 
or "two-dimensional and three-component" turbulence. Favier et a/.— also applied eddy- 
damped quasi-normal Markovian (EDQNM) model to study anisotropic turbulence in the 
quasi-static limit, and observed consistency with the numerical results. Burattini et a/.— 
showed that the nonlinear energy transfer at small scales is both radial and angular, thus 
necessitates an anisotropic treatment of the the quasi-static MHD turbulence. 

In this paper we derive an expression for the energy flux that includes the dissipation 
induced by the Lorentz force. We show that the energy flux has a strong wavenumber 
dependence that leads to the steepening of the energy spectrum. It is important to emphasize 
that we use isotropic Kolmogorov-like energy spectrum in the present model. The angular 
dependence of the spectrum would be incorporated in future. 

The organization of the paper is as follows. We describe the model and its results in 
Sections 2 and 3 respectively. Section 4 contains conclusions. 

II. THE MODEL 

The equation of motion for the liquid metal flow under the quasi-static approximation 

_ + „. Vu= _V Wp) __.__ + „VV (2) 

along with the incompressibility constraint V ■ u = (i.e., p is a constant). In the above 
equation, u is the velocity field, p is the pressure field, B is the mean magnetic field along 
the z direction, a is the electrical conductivity of the fluid, and v is the kinematic viscosity 



of the fluid. It is convenient to describe the flow properties in the Fourier space, for which 
we rewrite the Navier- Stokes equation in Fourier space as 

_ jL — = -ikip(k) J p - ikj^2uj(q)Ui(k - q) -(cos 2 9)ui - uk 2 u h (3) 

where 9 is the angle between the mean magnetic field and the wavenumber k. The interaction 
parameter N, which is the ratio of the Lorentz force and the nonlinear term, is aB^L/^pU). 
It is also important to keep in mind that the Lorentz force in the polar region is stronger 
than that in the equatorial region due to the cos 2 9 factor. Hence, the interaction parameter 
is an averaged measure. 

The corresponding energy equation is 

^!M = T (k) - 2^ cos 2 (#)£(k) - 2^ 2 £(k), (4) 

at p 

where E(\s) = |u(k)| 2 /2 is the energy spectrum, and T(k) is the kinetic energy transfer rate. 
The second and the third terms in the RHS are the dissipation rates due to the Lorentz force 
and the viscous force respectively. For N = and large Reynolds number Re = UL/u, the 
flow becomes turbulent, and the energy spectrum exhibits the famous Kolmogorov's k~ 5 ' 3 
power law in the inertial range. However for finite N, the Lorentz force induces additional 
dissipation that leads to the modification of the energy spectrum. The Lorentz force is 
dissipative in the quasistatic limit, so we can define another Reynolds number Re CT as the 
ratio of the nonlinear term pu ■ Vu and the Lorentz force (of dissipative nature), i.e., 

We term Re CT as the "Reynolds number based on resistivity" . We will show below that Re CT 
plays an important role in the dynamics of liquid metals MHD. 

Since the mean magnetic field induces anisotropy&i^£, the local energy density is best 
described using ring spectrum E(k, 9)^^, where k is the wavenumber of the ring, and 9 is 
the angle between the mean magnetic field and the wavenumber k of the ring, as shown in 
Fig. [TJ We extrapolate Pope's prescription^ for the energy spectrum in isotropic turbulence 
to the ring spectrum as 

E(k, 9) = C{Yi{k)) 2 "k-^f L {kL)Uk V ) 9 -^ (6) 
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FIG. 1. (a) Shell decomposition of the wavenumber space; (b) Ring decomposition in which a shell 
is further divided into rings. Adapted from Teaca et a/.—. 



where C ~ 1.5 is the Kolmogorov's constant, H(k) is the energy flux emanating from the 
wavenumber sphere of radius k, and g(9) represents the anisotropic component of the en- 
ergy spectrum. The functions fi,{kL)J n {krj) specify the components of the large-scale and 
dissipative-scale spectra respectively, and they have been prescribed by Pope^ as 



kL 



fL{kL) -\[{kLy + CL yi\ 
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(7) 
(8) 



where the cl,c v ,Pq, (3 are constants, which are determined by matching the above function 
with experimental results on isotropic turbulence at high Reynolds number. We take Cl ~ 
6.78, c v ~ 0.40, j3 ~ 5.2 and po = 2 as suggested by Pope^. We remark that the splitting 
of E(k, 9) into a product of functions of k and 9 is an assumption that needs to be verified 
by simulations and/or experiments. The quasi-static MHD turbulence is anisotropic for 
moderate to high JV's, hence the form of fi{kL) and f v (kr)) described by Eqs. (J7J) and (jSJ) 
may not be applicable to quasi-static MHD turbulence. Yet, in the absence of extensive 
data sets for such flows, we choose Pope's prescription for our calculations. In the present 
paper we focus on the inertial and dissipative range, fh{kL) = 1. Under the special case of 
isotropic turbulence, g{9) = 1, and 



d9E(k,9) = E{k) = C(n(A;)) 2/3 A;- 5/3 /r,(^) 



(9) 



In Kolmogorov's phenomenology of fluid turbulence, the energy flux Tl(k) is a constant in 
the inertial range since the dissipative term 2uk 2 E(k, 9) is negligible. In quasi-static MHD, 
the Lorentz force introduces an additional dissipation term. In contrast to the viscous 



dissipation, the dissipation due to the Lorentz force is active at all wavenumbers. As a 
consequence, the inertial-range energy flux II (A;) decreases significantly with the increase 
of k. Quantitively, the difference between energy fluxes Ii(k + dk) and Ii(k) is due to the 
energy dissipation in the shell (k, k + dk), i.e., 



or 



n(fc + dk) - n(fc) 
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- [ g{6)d9- c 2 = - [ g(6) cos 2 9d9. 
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(12) 



We integrate the above equation from k — k\, which is the starting wavenumber of the 
inertial range. Assuming that the energy flux at this wavenumber is Ho, we obtain 



n(fc) 



1/3 



l- 



2Cci 
2cic 3 C 



1/3 



h{krf) ~ 
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3p n 



2/3 
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h{kri) 



Hkv)-l^=Ni 2 (k v ), 



3 3 vc^Re 

where rj is the Kolmogorov length, and (i^/Ho?? 4 ) 1 / 3 = C3, or 

,3 \ 1/4 



(13) 



V 



-3/4 



(14) 



We choose C3 = 3.1 in order to achieve H(k) —> for fci]>l when N = (isotropic case). 
We also take Kolmogorov's constant C to be 1.5 for all our calculations. The constants c\ 
and c 2 depend quite crucially on g(9), which is a function of N. Unfortunately, at present, 
g{9) is not known accurately. Therefore, we choose c\ — 1 and c 2 = 1/2, which are the 
values for the isotropic or iV = case (g{9) = 1). The integrals ii and I 2 in the above 
equations are 

'•A;// 



h{kri) 
h(kv) 



kirj 
kij 

kiii 



dk'k ,1/3 f v (k') 

dk'k'- 5/3 f v (k'). 



(15) 
(16) 



Equation (TT3"j) indicates that the second term of the energy flux which arises from the Lorentz 



force depends crucially on the constant N' = 2c 2 CN/(3y/csRe), termed as the "normalized 
interaction parameter" . 
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Following Popei^, we assume that the inertial range wavenumber starts at around k\ = 
6 x 2%/L. Therefore, the lower limit of the integral is kit] = 6(2ir)(r)/L) = Ylix x (c3_Re) -3 / 4 . 
Note that the energy flux Tl(k) peaks at k = k\ with value n . In the limit v — > 0, and for 
fc's in the inertial range, we can approximate the energy flux as 

"n(&)l 1/3 _ „ c 2 CN 



n n 



[(k lV )- 2 / 3 - (k V )-^} (17) 



\Jc3Re 

which demonstrates that H(k) decrease as (A — Bk~ 2 ^ 3 ) 3 in the inertial range. Wavenumber- 
dependence of the energy flux is consistent with the numerical findings of Ishida and 
Kaneda 1 ^. 

The energy spectrum can be easily computed from the aforementioned fc-dependent en- 
ergy flux using 



curk^ 3 f v (k v ) 



n(fc) 



if h ^> hi 

E(k)=--{ ' L °° J (18) 

Clf /3 k^ 3 f L (kL), otherwise. 

Given the energy spectrum, the dissipation rates due to the viscous and Lorentz forces could 
also be computed as D v = J 2vk 2 c\E{k)dk and D^ = J (2aB 2 /p)c 2 E(k)dk respectively. 

III. RESULTS 

After the above discussion on the variable energy flux and the energy spectrum, we com- 
pute these quantities for two sets of parameters: Re = 10 4 and N = 0,1, 5; and Re = 2500 
and N = 0,1,5. We test whether our model provides results consistent with earlier ex- 
periments and numerical simulations. For R = 10 4 and iV = 0,1,5, Figure [2(a) exhibits 
the normalized energy spectra E(krj){kri) a / Eikirj), where a is the spectral exponent. Fig- 
ure [2(b) shows the corresponding normalized energy fluxes n(/c)/Ho (solid curves) and the 
approximate energy fluxes (dashed curves) computed using Eq. (IT7j) . These figures show 
that the energy spectrum follows a power law in the inertial range. Also, the approximate 
energy flux computed using Eq. (|T7|) matches very well with the actual energy flux computed 
using Eq. (1131) for kr\ in the inertial range. 

Figure [3] exhibits corresponding the energy spectra and fluxes for Re = 2500 and N = 
0, 1, 5. As evident from the normalized energy spectra plots, power law fits reasonably well 
with the energy spectrum for a relatively narrower inertial range. The approximate energy 
flux formula [Eq. (|T7|) ] do not match with the energy flux of Eq. (Tl3l) because the flows with 




FIG. 2. (a) Plots of normalized energy spectra E(krj)(krj) a / E{kivj) for Re = 10 4 and N = 
(red thin curve), 1 (green thicker curve), and 5 (thickest blue line). Here a = 5/3, 1.8, and 2.1 are 
the best fit spectral exponents, (b) Plots of the corresponding energy fluxes (solid lines) and the 
approximate energy fluxes computed using Eq. (|17p (dashed lines). 



Re = 2500 have a relatively smaller inertial range. Also note that the magnitude of the 
spectral exponents a for Re = 2500 is higher than that for Re = 10 4 because of the increase 
in the dissipation for lower Re. 



We compile the spectral indices and the ratio of the dissipation rates D^ / D v for all our 
runs in Table 1. The dissipation rate due to the Lorentz force (-Dv) is significant, in fact, it 
is more than D v for the parameters studied in the paper. We should however keep in mind 
that the system is anisotropic, and the angular dependence is important. For N > 1, the 
dissipative Lorentz force is stronger than the nonlinear term, or Re CT < 1. Yet, E(k) exhibits 
power law scaling possibly because the Lorentz force is active at all scale. The spectral 
index, as well as the the ratio of the dissipation rates Dn/D u , grow monotonically with the 
interaction parameter N. This trend holds qualitatively in experiments and simulations. 




FIG. 3. (a) Plots of normalized energy spectra E{krf){krf) a / E{kirf) for Re = 2500 and N = 
(red thin curve), 1 (green thicker curve), and 5 (thickest blue line). Here a = 5/3, 1.9, and 2.2 are 
the best fit spectral exponents, (b) Plots of the corresponding energy fluxes (solid lines) and the 
approximate energy fluxes computed using Eq. (|17p (dashed lines). 

IV. CONCLUSIONS AND DISCUSSIONS 



In this paper we construct a model based on variable energy flux for modeling liquid- 
metal MHD under quasi-static approximation. The energy flux decreases with the increase 
of wavenumbers due to the dissipation induced by the Lorentz force. Since this dissipation is 
effective at all scales, the energy spectrum in the inertial range gets steepened significantly. 
The magnitude of the spectral exponent increases monotonically with the interaction pa- 
rameter N. Our spectral exponents are in qualitative agreement with the numerical*^^ 
and experimental results^. 

Based on experiments and numerical simulations, some researchers argue that the energy 
spectrum of the quasi-static MHD is proportional to k~ 3 , similar to the forward enstrophy 
cascade regime of two-dimensional fluid turbulence. However, Favier et alM- showed that the 
velocity profile of quasi-static MHD turbulence differs significantly from 2D fluid turbulence, 
and it resembles "two-and-half-dimensional-" , or "two-dimension three-component (2D-3C)" 



TABLE I. Table depicting the spectral exponent a and D^/D^, the ratio of dissipation due to the 
Lorentz force and the viscous force, for various Reynolds number Re and interaction parameter 
N. The normalized interaction parameter N' = (2/3)(c2CN)/\/c3Re, where C is the Kolmogrov's 
constant, Re is the Reynolds number, and C2 is a parameter. Re CT = 1/N is the "Reynolds number 
based on resistivity" . 



Re N N' Re a a D N /D V 

10000 - 5.0/3 

10000 1 0.0028 1 1.8 1.8 

10000 5 0.014 0.2 2.1 15.0 

2500 - 5.0/3 

2500 1 0.0057 1 1.9 1.35 

2500 5 0.028 0.2 2.2 9.7 

turbulence with all the three components of the velocity field being important. We remark 
however that even for the aforementioned 2D-3C turbulence, the spectral exponent could 
be greater than 3 due to the additional dissipation induced by the Lorentz force, similar 
to the steepening of the energy spectrum for the two-dimensional turbulence with Ekman 
frictior*^^. Note that for fluid turbulence, the energy is conserved in the nonlinear triad 
interactions, and the dissipation takes place at small scales only. Thus, the quasi-static 
MHD turbulence in which dissipation takes at all scales is very different from pure fluid 
turbulence. A simple extrapolation of pure fluid turbulence to quasi-static MHD turbulence 
is not appropriate. 

The system under investigation is anisotropic due to the cos 2 (#) term in the Lorentz 
force term. We believe that the g(6) term of Eq. ([6]) could capture the anisotropic effects 
in a significant manner. Unfortunately our present calculation assumes isotropic energy 
spectrum (g(6) — 1, or c\ — 1 and C2 = 1/2) due to an inadequate knowledge of g(0). We 
are in the process of computing g(6) using direct numerical simulations, which will enable 
us to compute precise values of c\^ for various values of the interaction parameters. We 
expect the spectral exponents computed by this modified model to agree better with the 
numerical simulations and experiments. For time being, the isotropic energy spectrum is 
expected to be applicable for moderate iV's (N ~ 1). We need to extend the present analysis 
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to anisotropic situations, especially for large interaction parameters. 
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